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ABSTRACT 

A three dimensional finite element code is written for 
truss analysis and design. Trusses may be designed for mini- 
mum weight subject to constraints on: member stresses, Euler 

buckling, joint displacements and system natural frequency. 

The optimum configuration may be found in addition to optimi- 
zation with respect to member sizes. 

The finite element code may be used as a stand alone 
analysis tool or may be coupled to an optimizer of the user’s 
choice. The finite element displacement method of analysis 
is used for static analysis and eigenvalues are calculated 
using the subspace iteration technique. 

A major portion of this work is the software (user guide) 
which is presented in detail with examples and results . Ex- 
planation of how this code is coupled to an optimizer is 
given . 

Design variables are member area sizes, joint coordinates, 
or both. Examples are presented to demonstrate the method. 
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I . INTRODUCTION 



Design optimization requires the minimization or maximi- 
zation of some parameter. Optimization of structures has had 
continuing changes since its development in the early 1960's 
with an active area of research being elastic truss structures. 
The main goal is to design structural systems that efficiently 
perform specified purposes. For any design to be acceptable, 
it must satisfy a variety of physical, aesthetic and economic 
constraints . 

Since most physical problems can be modeled by some mathe- 
matical formulation a computer program can be written to per- 
form the necessary calculations . 

The purpose of this research was to develop a FINITE 
ELEMENT code that could be easily coupled to an optimizer, 
thus (1) allowing testing and comparison of various optimizers 
and (2) provide a useful design program in its own right. 

The design problem considered in this study is the opti- 
mization of three-dimensional indeterminate trusses , for 
multiple static load conditions. The objective is to mini- 
mize the weight of the structure where the design variables 
are member sizes and joint coordinates. Constraints include 
stress, displacement, Euler buckling, and the natural fre- 
quency of the system. 
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This document describes the use and capabilities of the 
finite element computer code to be coupled to an optimizer. 

The user's manual presented in Chapter V contains a simple 
design example in which the program is coupled to the CONMIN 
optimization code [Ref. ll. 

Additionally, guidelines for coupling the code to an 
optimizer of the user's choice are presented. 

Several examples demonstrating the program under a variety 
of conditions are presented. Conclusions and recommendations 
for future work are given. 
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II. OPTIMIZATION 



A. INTRODUCTION 

The main goal of structural engineering optimization is 
to design structural systems that efficiently perform speci- 
fied purposes. Selection of a specific algorithm must include 
the following considerations: 1) the structure should be 

analyzed as few times as possible, 2) the algorithm should 
minimize specific gradient information required, and 3) it 
should provide reasonable assurance that an optimum design 
will be reached. 

The next few sections discuss the general formulation of 
the algorithm with respect to the above requirements . 

5. FORMULATION 

Minimum weight design of trusses is presented in the 
general form of a mathematical programming problem as follows: 

Minimize F(X) (Eq. 1) 

Subject to: 

gj (x) <_ 0 j=l,m (Eq. 2) 

x^ <_ <_ xV i = l,n (Eq. 3) 

where F(X) is the objective function, in this case, weight of 
the structure be minimized. X is the vector of independent 
design variables and here contains the member cross-sectional 
areas as well as the coordinates of the joints . The inequality 



12 



constraints, g^(x) < 0,j=l,m, must be satisfied for the design 

to be accepted as feasible. These include limits on stress, 

Euler buckling, joint displacements and the first fundamental 

• . o n 

frequency of the structure . Side constraints x^ and x^ are 
lower and upper bounds on the design variable. These may be 
treated as general inequality constraints . 



C. PARAMETERS 



The design variables, constraints, and objective function 
considered in the design process are discussed here. 

1 . Design Variables 

The design variables the user has available are member 
cross-sectional area, reciprocal of member cross-sectional 
area, joint coordinates, or both member areas (or their reci- 
procals) and coordinates. The cross-sectional areas are A^ , 
k=l,NE, where NE is the number of elements in the truss 
structure. The joint coordinates are , i=l,2,3 j=l, NJ, 
where i is the coordinate axis, j is the joint number and NJ 
is the total number of joints . Design variable linking is 
allowed in both member sizing and coordinate design variables. 

2 . Objective Function 

The objective function considered here is weight. 



Weight (W) 



ME 
= Z 
i = l 




(Eq. 4) 



where p is the material density, A is the cross-sectional area 
and is the length of the member. The truss may be made up 
of members of differing materials . 
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3 . Constraints 

i=l ,NE (Eq. 5) 

a . . 

Stress: 1 - 1 <_ 0 j=l,NLC 

a 

i=l ,NE (Eq. 6) 

a . . 

-ii - 1 <_ 0 j =1 ,NLC 

a 

where a is the stress in member i under load condition j . NE 
is the number of elements and NLC is the number of loading 
conditions . 

a ~ is the lower bound on the stress (maximum compressive 
stress) and a + is the upper bound on the stress. The upper 
and lower bounds on stress may be different for each member, 
but are taken to be the same for every loading condition. 

4 . Euler Buckling 

The stress at which Euler buckling occurs is given by 

o. E = KlAl 1 i=l,HE (Eq. 7) 

1 L . 2 

1 

where the subscript i corresponds to the member number, E^ 
is Young's modulus and is a constant depending on the 
cross-sectional geometry of the member. 

5 . Displacement 

Displacement limits are imposed at prescribed joints 
to create a constraint equation as follows : 
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(Eq. 8) 



U ijk 



^- 1<0 

u ijk 



(Eq. 9) 



+ 



where u^^ and are lower and upper bounds on the displace- 

ment at joint i in the coordinate direction j under loading 
condition k. 

6 . Frequency 

The first fundamental frequency of the structure is 
required to exceed the specified lower bound, so that 



1 - X/X <0 



(Eq. 10) 



Reference 
structural design 



2 is an excellent source for the basic 
formulation . 
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Ill . FINITE ELEMENT METHOD 



A. INTRODUCTION 

Several features are desirable when the finite element 
methods of analysis are used in design optimization. First 
the number of analyses for the structure should be kept to a 
minimum. Second, the amount of gradient information required 
during the design process should be reduced to shorten run 
times and computer storage requirements. Third, the user 
should be able to specify only that gradient information 
desired . 

B. ANALYSIS 

Initial problem formulation includes member sizing, mate- 
rial properties (which may be different for each member) , a 
set or sets of external loads, and specified support condition. 

The analysis for the stresses and deflections must satisfy 
the conditions of equilibrium of forces at the nodes and com- 
patibility of deformation. In this analysis the weight of the 
individual members are not part of the specified load conditions . 

Additionally for this analysis the following assumptions 
are made. Trusses will be treated as discrete elements, and 
each element will be treated as pin-connected with loads and 
reactions supported at the joints. 
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The general method of solution is as follows. The Dis- 
placement (Stiffness) method reference 2 considers the joint 
displacement components as the unknowns , and written in 
matrix notation for the general case is 

Ku = P (Eq. 1) 

where K is the global stiffness matrix, P is the vector or 
vectors of applied loads, and u is the vector of vectors of 
displacements . 

Once displacements at every node are known, the internal 
forces and stresses are calculated by applying the appropriate 
force-deflection relationships. 

When the system's natural frequency constraints are con- 
sidered the design process requires the solution of an eigen- 
problem. This solution will determine the natural frequencies 
and normal modes of the structure. For linear elastic struc- 
tures, the finite element approach leads to the following 
equation of motion with free vibration conditions, 

Mu + Ku = 0 (Eq. 2) 

where M is the global mass matrix, and ii is the linear accel- 
eration vector. This leads to an eigenvalue problem of the form 
IK- A ld=0, and is solved here by the subspace iteration method 
reference 2 . 
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C. GRADIENTS OF CONSTRAINTS 



Gradient computation of necessary functions in a design 
optimization process arise from the need for derivative infor- 
mation for efficient mathematical programming. For structural 
optimization the gradient of displacement with respect to the 
design variables, 3u/3x, is needed from which VG, is calculated 
for stress, displacement, and buckling. 

Consider the finite element method Ku = P . Taking the 
derivative of both sides with respect to the design variable 
, we have : 



3 



[Ku] 



(P) 



(Eq. 3) 



l 




Assuming the loads P are not a function of x^ , 



3x 



l 



(K) ] (u) + K — (u) = 0 
S ~ S ~ 



(Eq. 4) 



Finally, we arrive at 



3 



(u) = -K" 1 [— (K)] u (Eq. 5) 



3x 



l 



where K - ^ is the inverse of K. For efficiency K - ^" is not 



actually calculated. 
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IV. PROGRAM FEATURES 



A. INTRODUCTION 

Each computer code has its own special features with which 
the user should be familiar if the program is to be used 
effectively. The next few sections will discuss an overview 
of the code and a typical problem that might be solved using 
the program. This problem along with other numerical examples 
will be presented in detail with results in Chapter VI. 

The FINITE ELEMENT code was written to be used as a stand 
alone analysis program or an analysis code that could easily 
be coupled to an optimizer (of the user's choice) through 
simple modifications to the main driver program. 

With user supplied area and coordinates coupled with input 
control parameters, the analysis mode will calculate the weight 
of the truss structure. Design variables may be chosen as 
member areas, joint coordinates, or both. Additionally, 
gradient information will be calculated with respect to area 
variables, coordinate variables or both. Coupled to an opti- 
mizer the code will optimize the weight of the structure and 
print the final optimization information and gradient vectors 
of specified constraints. 

The following example of a 25-bar space tower presented 
in Tables (I-IV) shows some of the options the code contains. 
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Table I 



25 -BAR TOWER SAMPLE OUTPUT 



INPUT CONTROL PARAMETERS 



I NPUT PARAMETERS FOR STRUCTURAL ANALYSIS AND DESIGN ROUTINE, "SAD" 
25--3AR SPACE TRUSS (TEST CASE STR ESS /OI S ?./ 3UCKLI NG/FREU. ) 



CONTROL parameters 



NUMBER OF ELEMENTS, ME = 25 
TOTAL NUMBER CF JOINTS, NJ = 10 
JOINT CONSTRAINT VARIABLE, NCJ = 10 
NO, OF MATERIAL TYPES, NMT = 1 
NO. OF LOAD CCND., NLC = 2 
NO. CF EIGENVALUES, N6IG = 2 
NO. OF EIGENVALUES CALC., NEIGl= 2 
NO. OF FIXED MASSES, NF MASS = 2 
BUCKLING CONSTRAINT ID, NEUBC = 1 
NO. OF DISPL. CONST., NOSPLC = 4 
NO. OF FRFQ. CONST., NF RE Q = 1 
LUMPED MASS OPTION, LMASS = 0 
VAOB. CALC. CONTROL , IDVCLC = 2 



ACCELERATION DUE TO GRAVITY, GRAV = 0.38640E+03 

EIGENVALUE CONVERGENCE TOLORANCE, EPSEIG = 0. 100006-03 



JOINT 

JOINT 


COORDINATES 

X 


Y 


l 


I 


-0.37500E+02 


0.0 


0. 200 OOE+O 3 


2 


0.37500E+02 


0.0 


0 .200006+03 


3 


-0.375006+02 


0. 3 75QQE+02 


0 . 1 00 00E+03 


4 


0.375 00E + 02 


0.3 75 OOF + 02 


0. 1 0000E+03 


5 


0 • 375006+02 


-0. 37500F+02 


0.1 00 00E+03 


6 


-0.37500E+02 


-0. 3 75 00F+02 


0.1 QQOOE+03 


7 


-0. 100006+03 


0.1 00 0QE + 03 


0.0 


3 


0. 10000E+03 


0. 1 00 OOF +03 


0.0 


9 


0. 10000E + 06 


-0. 100006+03 


0.0 


10 


-0. 100006+03 


-0. 10000E+03 


0.0 



COORDINATE OESIGN VARIABLES 





DESIGN 


VAR I ABLE 




MULTIPLIER 




JOINT 


X 


Y Z 




X 


Y 


l 


1 


0 


0 


0 


0.0 


0.0 


0.0 


2 


0 


0 


0 


0.0 


0.0 


0.0 


3 


1 


2 


3 


-0. 1000E+0L 


0.10005+01 


0. LOOOE+OI 


4 


1 


2 


3 


0. IOC 06 + 01 


0.1000E+0L 


0. 10006+01 


5 


1 


2 


3 


0.1000E+01 


-0 • I 000E + G1 


0. 1000E+01 


6 


1 


2 


3 


-0. LOQOF + Ol 


-0. I000E+01 


0.1 Q00E+01 


7 


4 


5 


0 


-0.L000E+01 


0. 1000E + 01 


0.0 


3 


4 


5 


0 


0. LOOOE + Oi 


0.1000E+01 


0.0 


9 


4 


5 


0 


O.IOGOE+Ol 


-0. 1000E+01 


0.0 


10 


4 


5 


0 


-0 . 10 OOE+O 1 


-0.1000E+01 


0.0 



MATERIAL PROPERTIES 

CODE E RHO SIG-MIN S IG-MAX K-EULER 

l 0. 1000E+C8 0.1000E+00 -0.4000E+05 0.4000E+05 0.3927E+02 
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Table II 



2 5 -BAR TOWER SAMPLE OUTPUT 

JOINT DISPLACEMENTS CALCULATED BY SYSTEM ROUTINE "SAD" 



LOAD CONDITION 
D IS PLACEMENTS 



L 



JOINT 

1 

2 

3 

4 

5 

6 
7 
3 
9 

10 



1 



•0.2 19C5E— 02 
0.219 10E-02 
0.9 07 89 5—0 1 
0.9 1278E-01 
•0.907895-01 
■0.9L278E-01 
0.0 
0.0 
0.0 
0.0 



DEGREE OF FREEDOM 
2 

0.3801 7E+00 
-0.380175+00 
-0 . I 5964E-0 l 
0. L7510E-01 
0.15964E-01 
-0. 1751 IE-01 
0.0 
0.0 
0.0 
0.0 



-0.27099E-01 
-0.27099E-01 
-0.68752E— 0 1 
0.361 00E-G1 
-0.68752E-01 
0. 36100 E-01 
0.0 
0.0 
0 .0 
0.0 



LOAD CONDITION 
DISPLACEMENTS 



JOINT 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 



1 



0.201 26 E-01 
0.22911 E-0 1 
0. 99 53 IE-03 
0 .647335-02 
0.81 493E— 03 
0.66536E-02 
0.0 
0.0 
0.0 
0.0 



DEGREE OF FR c EDOM 
2 



0.38860E+00 
0. 388605+00 
0.2595 1E-G1 
0.2670 76-01 
0.244355-01 
0.25192E-01 
0.0 
0.0 
0.0 
0 .0 



- C. 21 023E- 01 
-0.32637 E-01 
-0.95653E-01 
-0. 10297E+00 
0.62374 E-01 
0.701 94E-01 
0.0 
0.0 
0.0 
0.0 



JOINT COORDINATES 



JOINT 


X 


Y 


L 


l 


-0. 37500E+02 


0. 0 


0.2000GE+03 


■> 


0 . 37500E + 02 


0. 0 


0. 20000E+03 


3 


-0 . 37500 E +02 


0.37500E+02 


0 .100006+03 


4 


0.375005+02 


0. 375005+02 


0. 10000E+03 


5 


0. 375005+02 


-0.375 00E + Q2 


0. 10000E+03 


6 


-0. 375005+02 


-0.3 7500E+02 


0.10000E+03 


7 


-0 . 10000E+03 


0. 100005+03 


0.0 


8 


0. 100005+03 


0. 100005+03 


0.0 


9 


0. 10000E + 03 


-0. 10000E+O3 


0.0 


10 

JOINT 


-0. 10000E + 03 
COORDINATES 


-0. 100 006+03 


0.0 


JOINT 


X 


Y 


L 


l 


-0 • 37500E+02 


0.0 


0.20000E+03 


2 


0 . 3 75 D0E + 02 


0. 0 


0.2 0000E+03 


3 


-0.3 7500E+02 


0. 3 75 ODE + 02 


C.10000E+03 


4 


0.375005+02 


0.3 75 005+02 


0 • 100 00E+03 


5 


0. 3 75 COE +02 


-0.375 00E+02 


0.10000E+03 


6 


-0 .375005+02 


-0.3 75 005+02 


0. 10000E+03 


7 


-0. 10000E+03 


0. 1 00 00E+03 


0.0 


8 


0. 10000E+D3 


0. 1 0000E + 03 


0.0 


9 


0 . 1 00 00 E +03 


-0. 100005+03 


0.0 


10 


-0. 10000E+03 


-0. 1 00 00E+03 


0.0 
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Table III 



25-BAR TOWER SAMPLE OUTPUT 



AREA/COORDINATE VARIABLES 



THERE ARE 8 

INITIAL VALUES 
0.20000E+01 



LOWER ROUNDS 
0.10000E-04 



UPPER BOUNDS 
0.50000E+02 



THERE ARE 5 

INITIAL VALUES 
0.375Q0E+02 



LOWER BOUNDS 
O.LOOOCE+OO 



UPPER 30UN0S 
0.50000E+03 



AREA DESIGN VARIABLES 

0.20000E+0 1 0.20000E+01 0. 20000E+01 0.20 OOOE+Ol 0 .20000E+01 

O.IOO0OE-O4 0.10000E-04 0.10000E-04 0.10000E-04 0.10000E-04 

0.50QQ0E+02 0 .5 0000E+02 0.5000GE+02 0.50000E+02 0.50000E+02 
COORDINATE VARIABLES 

0.37500E+02 0 . 10000E+ 03 0. 10000E+03 0. IOOOOE +03 
0.10000E+00 0.1 OOOOE+OO 0 . 10000E+00 O.IOOOOE+OO 
0. 50000E+0 3 0. 50000E+ 03 0.50000E+03 0.50000E+03 



JOINT DISPLACEMENT CONSTRAINTS 



DIRECTION 1=X, 2=Y, 3=Z , 0=RESULTANT 







LOAD 


LOWER 


UPPER 


NOOE 


OIR. 


COND. 


BOUND 


BOUND 


1 


1 


I 


-0.35006*00 


0.350CE+00 


2 


i 


l 


-0.35006*00 


0.35 OCE + OO 


1 


2 


I 


-0.35006*00 


O.35OOE + 0O 


2 


2 


1 


-0.3500E+00 


0 .35 OOE+OO 



FREQUENCY CONSTRAINTS 
LOWER BOUNDS IN CPS 



FREQUENCY 

NUMBER 

L 



LOWER 

BOUND 

0.1600E+02 
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Table IV 



2 5 -BAR TOWER SAMPLE OUTPUT 



FINAL OPTIMIZATION INFORMATION 



F INAL OPTIMIZATION INFORMATION 

OBJ * O* 267256c +03 

DECISION VARIABLES (X-VECTO?) 

I) 0* 14595E + 01 0. 742B5E + 00 0, 

7) 0.85953F+)0 0.6957LE+00 0, 

13) C. 6 82805 + 02 



CONSTRA 
L) 

7 ) 

13) 

L 9 ) 

25) 

31) 

37) 

43) 

6 9 ) 

55 ) 

61 ) 

67) 

73) 

79) 

85 ) 

9 i ) 

97) 
103) 

L 09 ) 
115) 

121 ) 
127) 

133 ) 
139) 
145) 

151 ) 
157) 



INT VALUES (G- 
-0.23866E-01 - 
-0. 292525-01 - 
-0. 10090E+01 - 
-0.896075*00 - 
— 0*9 L 9 3 2 5 * 30 - 
-0. 1 02 59E + 01 - 
-0.104915*01 - 
-0.754476*00 - 
-0*1 102 15*01 - 
-0.76690^*00 - 
-0.111455*01 - 
-0.1 0366F **01 - 
-0.105895*01 - 
-0.110185*01 - 
— 0. 923645+00 - 
-0.117265*01 - 
-0.806 7 7 E + 0 0 - 
-0.115185*01 - 
-0.78602600 - 
-C. 84000E+00 - 
-0.811206*00 - 
-0 • 1131 6E +01 - 
-C.11028EO1 - 
-0 . 122 34 F +0 1 - 
-0. 72096EO0 - 
-0.6 806 7E * JO - 
-0.118316*01 - 



VECTO P ) 

0. 993 I 75 + 00 
0. 29252F-01 
0.991905*00 
0. LIO J9F+01 
0.1 J8 07E*91 
0. 9 74 L 3E *00 
0. 950 8QE * 00 
0.124555*01 
0. 39 7 89F + 00 
0 • 1 23 3 IE + 01 
0. 88546E+00 
0. 963 435+00 
0.941115+00 
0.898 1 7E **00 
0. 107645+01 
0.3274GF+00 
0. 119 325 + 01 
0.848 15E+00 
0.121405+01 
0. 1 16 OOF +01 
0. I 138HE+01 
0.86b 40E+00 
0.897 218 + 00 
0. 77661E+00 
0. L279GE+01 
0.131936+01 
-0.81690E + 00 



-0 

-o, 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 



98392E+00 

39367E+02 



, L0068E + 01 
19707E+01 
► 100 35 E +01 
,63131 E +00 
, 8 69 05 E + 00 
* 1 0420 E +01 
, 10797E+01 
,849 206 + 00 
► 1 06 2 7 E + 01 
, b56 83 E + 00 
, 107035+01 
,103156+01 
, 10507E + 01 
, 10570E + 01 
, 95 7 23 E+ 00 
, 13066E+01 
, 96989 F. - 01 
, 1 70 9o E + 01 
, 17345E-04 
.383256+00 
.272 19 E + 00 
. 150 73 E + 01 
, 13962E+01 
. 168 47 E+0 1 
. L 44 73 E + 00 
,212376-01 
.156126+01 



THERE APE 7 ACTIVE CONSTRAINTS 

CONSTRAINT NUMBERS ARE 

1 7 8 111 120 126 153 

THERE ARE 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

TERMINATION CRITERION 
ITFR EQUALS ITMAX 

NUMBER OF ITERATIONS = 20 

OBJECTIVE FUNCTION WAS EVALUATED 

CONSTRAINT FUNCTIONS WERE EVALUATED 

GRADIENT CF OBJECTIVE WAS CALCULATED 

GRADIENTS OF CONSTRAINTS WERE CALCULATED 



60 

60 

20 

20 



0. 1 0425 E +01 
0. 46 9346 +02 



-0.10068E+0L 
-0.10160 E+01 
-0 • 7 1 5E +00 
-0 . 1 2 349 E +0 1 
-0. 12349E + 01 
-0 .714656+00 
-0.12113 £ +01 
— 0. o5 536E +00 
-0.655366+00 
-0. 12113E+01 
-0.1047BE+01 
-0. 10479E+01 
-0. 10 l 70E + 01 
-0. 101 70E+01 
-0.11 JL0E+01 
-0.35724E+00 
-0 *85 724 F +00 
-0. 11 OLOE+Cl 
-0. LI 95 IE +01 
-0. 74 8 77 E +00 
-0. 74877E + 00 
-0. 11951 E +01 
-0. 994166+00 
-0.909206+00 
-0 • 9941 6E+00 
-0. 909196 + 00 



TIMES 

TIMES 

TIMES 

TIMES 



0. 1 1 271E +01 
0.13895E+03 



-0.99 316E+00 - 
—0.98405 E +00 - 

— 0. 12 854E + 01 - 
-0.765076+00 - 
-0. 76507E+00 - 
-0. 1 2354E+01 - 
-0. 78870E+00 - 
-0. L3v46F+0l - 
-0.134466+01 - 
-0. 78 870E+00 - 
-0.952L6E+00 - 
-0. 952166+00 - 
-0.983 04 E+00 - 
-0.983 04 E+00 - 
-0 .899026 +00 - 
-0. 114285+01 - 
-0.1L428E+01 - 
— 0. 899026 *00 - 
- 0. 8 0493 E + 00 - 
-0.125 12E+01 - 
-0.12512E+01 - 

— 0. 8 0493E + 00 - 
-0. 10053E+01 - 

— 0. 10 908E + 01 - 
-0.10058E+01 - 
-0. 1 0 908 E + 01 - 



0.735675+00 
0. 736 36 E +02 



0. 19707E+Q1 
0.100636*01 
0. 536845*00 
0. 13813E+01 
0. 138L3E+01 
0. 5 3b 84c + 00 
0.11298E+01 
0. 7o8 3 2E + 00 
0. 7883 2E +00 
0 .112986+01 
0. 10412E+01 
0.10412E+01 
0.10095E+01 
0. 10095E+01 
0.1471 9E+0 1 
0. 332845+00 
0 • 3 3284E + 00 
0. 14719E+01 
0. 17520E + 01 
0.316506-01 
0.315516-01 
0. I 752GE + 01 
0 • 982Q9E+00 
0. 72168E+00 
0. 96209E+00 
0. 7216 8 E + 00 
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V. USER GUIDE 



A. INTRODUCTION 

In developing any computer code for engineering analysis, 
it is necessary to additionally develop concise, easily under- 
stood software. This USER GUIDE is writen to be easily 
followed assuming minimal FORTRAN knowledge . The format 
follows that of the optimization code, COPES /CONMIN , 
reference 1. 

This chapter is devoted to acquainting the user with the 
code and necessary input data. A simple 3-bar truss analysis 
is used as the example . 

B . GENERAL FORMULATION 

Each data card is set up to show the input data necessary 
with the 3-bar truss example underneath. Formats are of "IS" 
and "F10.0" type. "I" formats must be right justified, and 
"F" formats must have the decimal point. The number of cards 
read per data block is listed at the bottom of each block. 

C. CONSTRAINTS 

Constraints are calculated and stored in the G vector as 
listed in the following chart. The total number of constraints 
NCON = 2*NDSPLC + NFREQ + NE*NLC. When any of the constraints 
are missing from the G vector, all constraints are moved up. 
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For example, if there is no frequency constraint, then a dis- 
placement constraint would fill the first location of the G 
vector . 



r 




■\ 

i 



NFREQ 

(FREQUENCY CONSTRAINTS) 





1 



1 2* NDSPLC * NLC 

( DIS PLACEMENT CONSTRAINTS ) 










1 



1 



MEMBER #1 
LC # 1 



1 



1 



MEMBER # 2 
LC #1 



1 

J 



MEMBER #NE 
LC #NLC 



STRESS AND 
BUCKLING 



STRESS AND 
BUCKLING 

STRESS AND 
BUCKLING 
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D. EXAMPLE 



The initial layout of the 3-bar truss is shown in Figure 
6.1. Stress constraints were imposed as well as constraints 
on Euler buckling, displacement, and first fundamental 
frequency . 

With geometry specified as per the figure, 2 independent 
load conditions (P1,P2) were imposed with 3 member sizing 
variables (Al, A2 , A3) linked so that A1=X(1), A2=X(2), 

A3=X( 1) . 

1 . Proper ties /Conditions 

3 

The material used had a density (p)=0.1 lb/in . 

7 

Young's Modulus was selected as 10 psx. 

LOADS: P1 = P2 = 20,000 lbs. 

STRESS LIMITS: -15,000 < < 20,000 psi i=l,3 

Symmetry was maintained by linking of variables so 
A1=A3. Initial design began with A1=A2=A3= 1.0 in. One non- 
structural mass of 500 lbs was attached at joint 2 . Displace- 
ment constraints were imposed between -.02 and .02 inches in 
the x-direction with the first natural frequency limited to 

a value w < 20 Hz. 
n 

2 . Input Control Parameters 

The following input control parameters are given for 
ease of the following example. NE=3, NJ = 4, NCJ=4, NMT=1, NLC=2, 
NEIG = 2, NEIG1= 2 , NFMASS = 1, NEUBC = 1, NDSPLC = 2, NFREQ = 1 , LMASS = 0, 
IDVCLC=2. Table V is a listing of commonly used nomenclature. 
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TABLE V 



COMMON VARIABLE NOMENCLATURE 



A- MEMBER'S CROSS-SECTIONAL AREA 
BL- LOWER BOUND ON DISPLACEMENTS 

BO- OPPER BOUND ON DISPLACEMENTS 

DIR- DIRECTION 1=X, 2=1. 3=Z, 

0=LIM IT ON RESULTANT DISPL. 

E- YOUNGS MODULUS 
ELNO- ELEMENT NUMBER 

EPSEIG- CONVERGENCE TOLERANCE OF EIGENVALUE 
SOLUTIONS (D EFAULT = . 000 1) 

FC 1 , FC2 FCN- LOWER BOUND ON FIRS T. SECOND , ETC. 

SYSTEM NATURAL FREQUENCY 
FX- LOAD FORCES BEING APPLIED IN THE X DIRECTION 

FY- LOAD FORCES BEING APPLIED IN THE Y DIRECTION 

FZ- LOAD FORCES BEING APPLIED IN THE Z DIRECTION 

GRAV- ACCELERATION DUE TO GRAVITY (DEFAULT= 386. 4) 
IX- CONSTRAINT IDENTIFIER. IF NON-ZERO THE X-DOF 
IS CONSTRAINED 

IY- CONSTRAINT IDENTIFIER. IF NON-ZERO THE Y-DOF 
IS CONSTRAINED 

IZ- CONSTRAINT IDENTIFIER. IF NON-ZERO THE Z-DOF 
IS CONSTRAINED 
JN- JOINT NUMBER (GLOBAL) 

KEULER- EULER BUCKLING COEFFICIENT FOR BAR 
ELEMENTS 

LC- LOAD CONDITION 

LM ASS- LUMPED MASS OPTIONS. IF LMASS . NE. 0 LUMPED 
MASS MATRIX USED 

IF LMASS =0 CONSISTENT MASS MATRIX USED 
NEIG1-NUMBER OF EIGENVALUE/EIGENVECTOR TO BE 
EVALUATED DEFAULT = MIN(2*NEIG, NEIG+8) 

NEUBC- BUCKLING CONSTRAINT IDENTIFIERS . IF NEUBC. 
NE.O EULER BUCKLING WILL BE IMPOSED IN 
EAR ELEMENTS 

NFMASS-NUMBER OF FIXED MASSES ATTACHED TO THE 
STR UCTO RE 

NFREQ-NUMBER OF FREQUENCY CONSTRAINTS 

NID- NUMBER OF INDEPENDENT DEGREES OF FREEDOM 

NJ-NUMBER CF JOINTS 

NLC- NUMBER OF LOADING CONDTIONS 

NLJ- NUMBER OF LOADED JOINTS FOR THIS LOAD CONDITION 

NMT- NUMBER OF SEPARATE MATERIAL TYPES 

RHO- MATERIAL DENSITY 

SIGMIN- MINIMUM ALLOWABLE STRESS 

SIGMAX- MAXIMUM ALLOWABLE STRESS 

XA- INITIAL VALUE OF AREA DESIGN VARIABLE 

XAL- LOWER BOUND 

XA U- UPPER BOUND 

XC- INITIAL VALUE OF COORDINATE DESIGN VARIABLE 
XCL- LOWER BOUNDS 
XCU- UPPER BOUNDS 
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The following user's manual is divided into blocks A 
through M. Appearing directly below each data field line 
are the parameters for the 3-bar truss example . It is impor- 
tant to note the user may choose any units; however, all units 
must remain consistent throughout the problem. 
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r 



DATA BL OCK A 
DE SCRIPTIO N: Title Card 

Format and Example 




FIELD 



CONTENTS 



1 ANY 80 CHARACTER TITLE 

MAY BE GIVEN ON THIS LINE 



j 
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1 



DATA BLOCK B 

DESCR IP TIO N: Control Parameters 



Format and Example 



NEj NJ i 


NCJ 


NMT 


NLC 


NEOBCl NDSPLC 

| 


NPR I 


FORMAT 


1 

1 1 
1 

i i i 






i — j 


1 

1 

1 

L 1 1 


j 


815 

i— i 



j 3 ' ' 

1 1 1 J 


NJ 

to 

o 


1 






NEIG 

2 

J 


NEIG 1 
2 

1 J 


NFMASS 

1 


NFREQ 

1 

j 


I DVCLC| IRECP 

2 i 0 

1 

1 l_ 


FORMAT 

515 

1 




i 1 

1 

L , 


1 1 i r 1 

i iii 

i i i i i 

l ... . L, ... 1 1 1 1 



NOTE: DEFINITIONS OF INPOT CONTROL PARAMETERS FOR 
PROGRAM NEXT PAGE. 






j 



■ * 

FIE LD CONTMI 

1 NE-number of elements 

2 NJ-number of joints 

3 NCJ-number of constrained joints 

4 NMT-number of seperate material types 

5 NLC-number of load conditions 

6 NEOBC- buckling constraint identifier 

(If NEUBC.NE.O -EULER buckling con- 
straints will be imposed on bar ele.) 

7 NDSPLC-number of displacement constraints 

8 NPRl-input print control (if NPRI.ne.O 

input information will not be printed) 

1 NEIG-number of precise eigenvalues to 

be evaluated 

2 NEIGI-number of eigenvalues to be evaluated 

DEF AUL T=mir. . of (2*NEIG , NNEIG + 8) 

3 NFMASS-number of fixed masses attached to 

structure 

4 NFREQ-number of frequency constraints 

5 IDVCLC-desig n variable control parameter 

If (I DVCLC.EQ. 1) NDV=NDVAR 1 

If (I DVCLC. EQ. 2) NDV=NDVAR 1 ♦ NDVAR2 

If (I DVCLC.EQ. 3) NDV=NDVAR2 

6 IRECP-recripocal variable identifier 

If (IRECP.GT.O) the X-vector contains 
the recripocal of area and VLB S VUB 
are changed accordingly 

i 
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DATA BLOCK C 

DES CRIP T I ON: Dynamic Analysis Information 
Format and Example 



LMASS 




GRAV 

1 

j 

1 

1 


EFSEIG 

1 , 

1 

l 


T FORMAT 

j IS. 

| 2F10.0 
1 


i 1 

1 o 
1 

L J 


386.4 | 
1 

L — — — ■ 1 


r o. i 

i i 

L_ J 


1 

i 

1 



FIELD 



CO NTEN TS 



1 LMASS-lumped mass option (if LMASS. NE.O) 

the lumped mass matrix is used. 

2 GRAV-accleration due to gravity 

(default =386 . 4 inches/sec ) 

3 EPSEIG-conv erge nee tolerance on eigenvalue 

solution, (def ault= . 000 1) 
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DATA block d 

DESCR IPTIO N; Joint Coordinates & 

D esign Variable Linking 

Format and Example 



55 

l — — — — — - 


X 


Y 


Z 


IX 


IY 


IZ 


PCX 


PC Y 


PCZ 


FORMAT 

15, 

3F10.0, 
315, 
3F10. 0 


L 1 


L— J 


1 


L 1 


1 1 


i j 


i i 


L — J 


L 1 


1 1 


—j 




1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



JN- joint coordinate number 
X-x coordinate 
Y-y coordinate 
Z-z coordinate 

IX-design variable associated 
IY-design variable associated 
IZ-design variable associated 
PCX-partici pati cn coefficient 
PCY-partici pa ti on coefficient 
PCZ-partici pati on coefficient 



with x coord, 
with y coord, 
with z coord, 
of x-coord. 
of y-coord. 
of z-coord. 



NOTE: Number of cards read=NJ 
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DATA BLOCK E 

DES CRIPTIO N; Material Properties 
Format and Example 





— — n 


r — i 


— r 1 


E 


RHO 


SIGMI N 


SIGMAX | K EULER | FORMAT 

j | 








I J 5F10.0 

1 1 


i i 


i , i 




L. — 1 i 1 



riOOOOOO. 

1 




. 1 


-15000, I 


20000.1 


4. 




L , 


I 




1 

L ^ 


1 

1- 




1 

1 



FIELD 

1 

2 

3 

4 

5 



NOTE: Number of cards read=NMT 

S 



CONTENTS 

E-Young's Modulus 
RHO- material density 
SIGM IN-minimum allowable stress 
SIGM AX-maximum allowable stress 
KECJL ER-Euler buckling coefficient 
for bar element 
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DATA BLOCK F 

D ES C RIPTION : Bar Element Information 
Format and Example 



ELNO | 


I node i I node 2 

i i 


| MATCOD 


NDSG 


1 

k 


FORMAT 




! 

! 


1 

i 

-J 






515, F 1 0 .0 J 


L 


_ 


1 


i i 



FIELD 



1 

2 



4 

5 



[ ~n n 

i i i 


p — 1 

i i 


1 

1 i 


i 1 

1 J 


H 7 ! 

1 1 




1 2 i 3 

I l 


2 

i... i 


1 I 1 

1 ! 2 ! 

1 L J 


1 '•! 
1 J 





4 j 2 




i i 
1 




’•j 


1 


1 

1 — L 


j 


_J 


J 


1 



C ONTE NTS 

ELNO-element number 
NODEI-global number associated with 
element node 

NODE2-global number associated with 
element node 

MATCOD -material type of this element 
NDSG-design variable number associated 
with this element 
A-member cross-sectional area 



NOTE: Number of cards read=NE 
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DATA BLOCK G 



D ESCRIPTION : Joint Constraint Information 



Format and Example 



i ix 

i 


IY 


IZ 


* 


FORMAT 










415 


] 








j 



' iii 

. J , 


i 

j 


1 

i i 


i 1 

1 

I 

i 




2 i 0 jo 

1 1 1 

1 1 1 1 


1 

i i 


1 

J 




1 3 1 
L J 


i 1 i 

i i 


L 1 


1 

1 


n 

i j 




i i 


i 

i 


1 


1 





FIELD CO NTEN TS 

1 JN- joint number 

2 IX x r y,z constraint identifier, (if non- 

3 IY zero the corresponding degree of freedom 

4 IZ is constrained) 

NOTE: Number of cards read=NCJ 
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DATA BLOCK H 

OMIT this card if NLC=0 was read in block B. 
DE SCRIPTIO N; Joint Loading Information 

Format and Example 



r nlj r 



| FORMAT I 



15 



r — t 

i 



i — r 



JN | FX 


FY 


FZ 


1 

1 






1— 1 


j 


. 



I FORMAT | 

I 15 , 3F 1 0 . 0 | 

I I 

1 



F 2 I 14140. | -14140.1 0~, 1 1 

II I I II 

I 1 1 J 1 I 

r 



2 


-14 140. 


-14140. 


0. 


1 


i i 


L 1 


L L 



FIELD 



CO NTEN T 



1 NLJ -number of loaded joints 

for this load condition 

1 JN- joint number 

2 FX 

3 FY- Forces in the X,Y,Z directions 

4 FZ 
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DATA BLOCK I 

OMIT this block if NFMASS=0 was read in block B. 
DE SCRIPTION : Lumped Mass Information 



Format and Example 



J25 

. — J 


MASS 
1 


1 

1 

J 


| FORMAT 

|I5,F10 . 0 
1 _ 1 






l ! , 


500. | 

1 

1 


1 



FIELD 



CONTENTS 



1 JN-joint number 

2 M ASS -concentrated mass at joinr (JN) 

mass is in force units 



_ , , ... i 
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DATA BLOCK J 

OMIT this block if NDVAR1=0 

DESCR I PTIO N: Design Variable Information 

(AREA Variables) 

Format and Example 





XA (I) ] 

| 


r r 


FORMAT 




| 


1 1 

i i 

1 L. 


8F10. 0 


L 


1 


» 



1 . 



1 . 



.01 

J 


. 01 

L 


1 

1 

I 


\ 

i 

i 


1 XAL (I) 
1 












10. ] 

1 


1 10. 

1 


1 

1 

L 


1 

i 

— j 


n xau (i) 

1 

J 



FIELD 



CONTENTS 



1 XA- initial value of area design variables 

2 XAL -lower bounds on area design variables 

3 XAU-upper bounds on area design variables 



NOTE: read one value of XA, XAL, XAU for each 

independent area variable defined in Block D 



Number of cards read =as required 
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DATA BLOCK K 

OMIT this block if NDVAR2=0 

DESCR IPTION: Design Variable Information 

(COORDINATE Variables) 



Format and Example 



XC (1) 



XC<2) 



r — ToTT 



r CoT 



XC (NDVAR1) 



r 20. 


r 

1 


1 . 


L 



FORMAT 
8F10. 0 



I XCL 



leu 



FIELD 



1 

2 

3 



CONTENTS 

XC-initial value of coord, design variables 
XCL-lower bounds on coord, design variables 
XCD-upper bounds on coord, design variables 



NOTE: read one value of XC, XCL, XCU for 

each independent coordinate variable 
defined in Block D 

Number of cards read =as required. 
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DATA BLOCK L 

OMIT this block if NDSPLC=0 in Block A 
D ES CR IPTIO N: Joint Displacement Constraint Information 

Format and Example 



JN 


DIR 


LC 


BL | BO 


| FORMAT | 

I 1 








1 

1 

I 


| 315, 2 F 10.0 


1 




. 


1 

L 


A - 1 



I 

L 

r 

L. 



2 I 

l 

2 r 

L 



i 



2 



1 


1 


-.02 


TT — 1 

.02 | 

1 




1 




j 


1 

L 1 






T -. 02]" .02 

1 L 



FIELD 



CONT ENTS 



1 JN-joint number 

2 DIR-direction 1=X 2=Y 3=z 

0=limit on resultant displ. 

3 LC-load condition 

4 BL -lower bound on displacement 

(If DIR=0 read 0 here) 

5 30-upper bound on displacement 



Number of cards read= NDSPLC 

, 
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DATA BLOCK M 



OMIT this block if N FREQ=0 was read in Block A 
DES CRIPTIO N: Frequency Constraint Information 

Format and Example 



FC 1 


FC2 


FC3 . 


f 1 

i 

j 


1 

I 


I 

1 


| FORMAT 
j 8F10. 0 




1 

1 

j 


20. | 


1 

! i 






1 i 

1 1 



FIELD 



CONT ENTS 



1 FC1- lower bound on first natural frequency 

constraint in CPS. (cycles per second) 
N FCN- lower bound on NFREQ-th natural 

frequency constraint in CPS. 

(cycles per second) 



NOTE: OMIT this block if NFREQ=0 was read in data 
block B. 

Number of cards read= as required 

1 
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VI. NUMERICAL EXAMPLES 



A. INTRODUCTION 

Design of planar trusses and 3-dimensional space towers 
are presented here and the corresponding numerical results 
are summarized to demonstrate the use of the code . 

In the examples given here, the design variables are 
member cross-sectional areas and joint coordinates. It 
should be emphasized that in practical design, the reciprocal 
of the member areas is usually a better choice for the design 
variables. However, the purpose here is to knowingly create 
difficult optimization problems, thus the choice of variables. 

The examples begin with the 3-bar truss . 

B. CASE 1: 3-BAR PLANAR TRUSS 

The simple 3-bar planar truss, as shown in Figure 6.1 has 
been previously used for the user guide example. This struc- 
ture was designed for optimum geometry subject to a set of 
two load conditions, buckling constraints, displacement con- 
straints, and a lower bound on the system's first natural 
frequency. The allowable stresses specified are 

-15000. < < 20000. psi i=l,3 

The member areas were linked in the following groups: A1=A3; 

A2 . 
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Case la 



1 . 

Only the stress constraint was imposed for this case. 
Final design information is given in Table VI. The number of 
analysis for this design was 38. 

2 . Case lb 

This case includes the previous constraint plus con- 
straints on buckling ( KEULER=4 . 0 ) , displacement (-.02 -.02 in.), 
and first natural frequency ( =20 . 0 ) . Results are given in 
Table VII. The number of analysis required for this design 
was 16 with 1.73 seconds of CPU time. 




Figure 6.1 3-BAR TRUSS 
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UU1 M 



C. CASE 2: 18-BAR PLANAR TRUSS 

A cantiliver truss 3 as shown in F igure 6 . 2 ? has been pre — 

viously used as a standard test case for structural design 

[Ref. 2 ] . The structure was analyzed for a single set of load 

conditions with allowable stresses being 

-20000. < a i < 20000. psi i=l,18 

1 

ioung's modulus was taken as 10 psi with a material density 
p = 0 . 1 lb . /cu in . 



250" 



250" 



250" 



250" 



250" 




Figure 6.2 18-BAR TRUSS 
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The independent coordinate variables were taken as X3, Y3, 

X5 , Y5, X7, Y7, X9 , Y9 . The member areas were linked as 
follows: A1=A4=A8=A12=A16 ; A2 =A6 =A10 =A14 =A18 ; A3=A7=A11=A15 ; 

A5=A9 =A13=A17 . There are a total of four independent area 
variables and eight coordinate variables . 

1 . Case 2a 

This design analysis includes stress, buckling (KEULER= 
4.0), displacement (-10.0 to 10.0 in.), and a first fundamental 
frequency of 3 Hz. Additionally, a nonstructural mass of W=5000 
lbs. was attached to node 1. The number of analyses for this 
design was 62. Results are presented in Table VIII. 26.72 
seconds of CPU time were used for this calculation. 

D. CASE 3: 2 5 -BAR SPACE TOWER 




Figure 6.3 25-BAR SPACE TOWER 
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The 25-bar space tower shown in Figure 6.3 was designed 
for two independent load conditions given in Table X. The 
allowable stresses were specified as follows 

-40000. < < 40000. psi i=l,25 

7 

Young’s modulus was selected as 10 psi with a material density 
p=.l lb/cu. in. Members are assumed tubular with a nominal 
diameter to thickness ratio of D/t=100 giving KEULER= 39 . 2 74 . 
Symmetry was imposed in both the x-z and y-z planes. Non- 
structural masses of W=500 lbs. were attached at nodes 1 and 
2. Coordinate variables were X4 , Y4 , Z4 , X8, and Y8 with the 
remaining coordinates linked to maintain symmetry. Area vari- 
ables were linked in the following manner: A1 ; A2=A3=A4=A5; 

A6=A7=A8=A9; A10=A11; A12=A13; A14 =A15 =A16 =A1 7 ; A1 8=A19 =A20 =A21 
A2 2 =A2 3=A2 4 =A2 5 . 

1 . Case 3a 

Stress, displacement (-0.35 to 0.35 in.), Euler buck- 
ling and first natural frequency limited to a value > 16 
Hz. were imposed for this case. Final design information is 
given in Table XI. The number of analyses for this design was 
61 with 27.99 seconds of CPU time used. 

E. CASE 4: 234-BAR SPACE TOWER 

The inital layout of the tower is shown in Figure 6.4 
stress limits were as follows 

-15000. < < 20000. psi i=l,234 

7 

Young's modulus was chosen as steel 3.xl0 psi with a material 
density of aluminum p =.l lbs/cu. in. 



1 . 



Case 4a 



Stress constraints as well as constraints on Euler 
buckling were imposed. The resulting problem had 56 area 
design variables and 42 coordinate variables. This problem 
required 91 structural analyses using 72 minutes of CPU time. 
The objective function was evaluated 89 times and gradients 
were calculated 30 times. Results are presented in Tables XII 
through XIV. Although the weight of the structure was only 
reduced by 2%, it should be noted that constraints were ini- 
tially violated and the optimizer overcame the constraint 
violations . It is believed the optimum that was reached is 
not the true optimum because of the extreme non-linearity of 
the problem. 
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FRONT/BACX 



© 



LEFT/RIGHT 



/ 

Figure 6 . 



4 2 34 -EAR SPACE TOWER 
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Table VI 



3-BAR TRUSS. DESIGN INFORMATION (STRESS) 



BODY 


LOAD 


JOINT 




LOADS 


COND. 




X 


Y Z 


1 


2 


0 . 1414E + 0 5 


-0.1414E+05 0.0 


2 


2 


-0 . 1414E+05 


-0.1414E+05 0.0 


FIXED MASSES 


NO . 


JOINT 


MASS 




1 


2 


500 lbs. 




MEMBER 




AREA (sq. in.) 
INITIAL 


FINAL 


A1=A3 




1.0 


. 7662E + 00 


A2 




1.0 


. 4 8 3 3E + 00 


JOINT 




COORDINATES (in.) 
INITIAL 


FINAL 


2 




10 


. 9 856E + 01 



FINAL WEIGHT FOR AREA AND COORDINATES 
WEIGHT = 2.6349 (lbs.) 
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Table VII 



3-BAR TRUSS. DESIGN INFORMATION 
(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 



BODY 



JOINT LOADING INFORMATION 
SAME AS TABLE VI 



JOINT DISPLACEMENT INFORMATION 



JOINT DIR. LOAD 

NUMBER COND. 



LOWER 

BOUND 



UPPER 

BOUND 



2 



- .02 



.02 



2 



1 



-.02 



.02 



MEMBER 
A1 = A3 
A2 



AREA (sq. in.) 
INITIAL 
1.0 
1.0 



FINAL 

0 . 99995E + 00 
0 .1000 E-01 



JOINT 

2 



COORDINATES (in.) 
INITIAL 
10 .0 



FINAL 

0 . 999 86E + 01 



FINAL WEIGHT = 2.8381 (lbs.) 
CRITICAL CONSTRAINT - FREQUENCY 
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Table VIII 



18-BAR TRUSS. LOAD CONDITIONS 

BODY 



LOAD CONDITIONS (lbs.) 



LOAD COND. 


JOINT 




LOADS 








X 


Y 


Z 


1 


1 


0.0 


-0 .2000E + 05 


0 .0 


1 


2 


0.0 


-0 .2000E + 05 


0 .0 


1 


4 


0.0 


-0 . 2000E + 05 


0 .0 


1 


6 


0 .0 


-0 . 2000E + 0 5 


0 .0 


1 


8 


0.0 


-0 .2000E+05 


0 .0 





FIXED MASS INFORMATION 




NO . 


JOINT 


MASS 


1 


1 


5000 (lbs 
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Table IX 



18-BAR TRUSS. DESIGN INFORMATION 
(STRESS, DISPLACEMENT, EULER BUCKLING, FREQUENCY) 



AREA ( sa . in . ) 



MEMBER 


INITIAL 


FINAL 


A1=A4=A8=A12=A16 


10 .0 


38.20 


A2=A6=A10=A14=A18 


10.0 


38.27 


A3=A7=A11=A15 


10 .0 


15.98 


A4=A9=A13=A17 


10 .0 


20 .62 




COORDINATES (in.) 




JOINT 


INITIAL 


FINAL 



X Y X Y 



3 


1000 . 


0 . 


859 . 


0 


5 


750 . 


0 . 


628 . 


0 


7 


500 . 


0 . 


384 . 


0 


g 


250 . 


0 . 


148 . 


0 



FINAL WEIGHT = 13,940 (lbs.) 
CRITICAL CONSTRAINT - FREQUENCY 
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Table X 



BODY 



2 5 -BAR TRUSS. LOAD CONDITIONS 



LOAD CONDITIONS (lbs.) 



LOAD COND. 


JOINT 


X 


LOADS 

Y 


Z 


1 


1 


0 .0 


2000 . 


-5000 




2 


0 .0 


-2000 . 


-5000 


2 


1 


1000 . 


10000 . 


-5000 




2 


0 . 


10000 . 


-5000 




3 


500 . 


0 . 


0 . 




6 


500 . 


0 . 


0 . 





FIXED MASS INFORMATION 




NO . 


JOINT 


MASS 


2 


1 


500 . (lbs . 



2 500. (lbs. 
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Table XI 



25-BAR TRUSS. DESIGN INFORMATION 
(STRESS, DISPLACEMENT, EULER BUCKLING, FREQUENCY) 



AREA (sq. in.) 



MEMBER 


INITIAL 




FINAL 




A1 


2.0 




1.46 




A2=A3=A4=A5 


2 .0 




0.74 




> 

CD 

II 

> 

II 

> 

00 

II 

> 

CD 


2 .0 




0.98 




A10 =A11 


2 .0 




1.04 




A12 =A13 


2 .0 




1.12 




A14=A15 =A16 =A1 7 


2 .0 




0.74 




A18=A19=A20=A21 


2 .0 




0 . 86 




A2 2 =A2 3 =A2 4=A2 5 


2.0 




0 .70 






COORDINATES (in.) 








JOINT 


INITIAL 




FINAL 






X Y Z 


X 


Y 


Z 


4 


37.5 37.5 100 


39 .4 


46 .9 


138.9 


8 


100.0 100.0 0.0 


73.6 


68 . 3 


0.0 



FINAL WEIGHT =267 (lbs.) 
CRITICAL CONSTRAINT - FREQUENCY 
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Table XII 



BODY 



234-BAR SPACE TOWER. LOADING INFORMATION 





LOAD 


CONDITIONS (lbs 


.) 


LOAD COND. 


JOINT 


X 


LOADS 

Y 


1 


49 


6000 . 


-20000 . 


1 


50 


6000 . 


-20000 . 


2 


49 


6000 . 


-20000 . 


2 


50 


-6000 . 


-20000 . 


3 


49 


6000 . 


-20000 . 


3 


50 


3000 . 


-10000 . 


4 


49 


3000 . 


-10000 . 


4 


50 


3000 . 


-10000 . 


5 


49 


-3000 . 


10000 . 


5 


50 


-3000 . 


10000 . 





FIXED MASS INFORMATION 




NO . 


JOINT 


MASS 


1 


49 


200 . 


2 


50 


200 . 



Z 

0.0 
0.0 
0 .0 
0.0 
0 .0 
5000 . 
-5000 . 
5000 . 
5000 . 
-5000 . 



(lbs . ) 
(lbs . ) 
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Table XIII 



234-BAR SPACE TOWER. DESIGN INFORMATION 



AREA (sq. in.) 



MEMBER INITIAL 

A1,A2,A3,A4 25.0 

A5,A6,A7,A8 25.0 

A9 , A10 , All , A12 25 .0 

Ak3 , A14 , A15 ,A16 25.0 

A17,A18,A19,A20 25.0 

A21,A22 ,A23,A24 25.0 

A2 5 , A26 , A2 7 , A2 8 25 .0 

A29 ,A30 ,A31,A32 25.0 

A33,A34 ,A35 ,A36 25 .0 

A37 ,A38 ,A39 ,A40 25.0 

A41 , A4 2 , A4 3 , A44 25.0 

A4 5 ,A46 ,A4 7 , A4 8 25 .0 

A49 ,A50 ,A51,A52 25.0 

A5 3 , A5 4 , A5 5 , A5 6 25 .0 

A5 7 ,A5 8 , A5 9 ,A60 25 .0 

A61,A62 ,A63,A64 25 .0 

A6 5 , A66 , A6 7 , A6 8 25 .0 

A69 ,A70 ,A71,A72 25.0 

A73,A74 ,A75 ,A76 25.0 

A77 ,A78 ,A79 ,A80 25.0 

A81 , A8 3 25 .0 

A82 , A84 25.0 

A85 ,A87 25.0 

A86 , A8 8 25 .0 

A89 , A9 0 , A9 1 , A9 2 , A9 3 , A9 4 ,A9 5 , A9 6 25 .0 

A9 7 ,A9 8 ,A9 9 ,A100 ,A101 ,A102 ,A10 3 ,A104 25 .0 

A150 ,A106 ,A10 7 ,A10 8 ,A109 ,A110 ,A111 ,A112 25 .0 

A113 ,A114 ,A115 ,A116 , A117 , A118 ,A119 ,A120 25 .0 

A121 , A122 , A12 3 ,A124 , A12 5 , A12 6 , A12 7 ,A12 8 25 .0 

A129 , A1 30 , A1 31 , A1 32 ,A133 ,A134 ,A135 ,A136 25.0 

A1 37 ,A138 ,A139 ,A140 , A141 ,A14 2 , A14 3 , A144 25.0 

A14 5 ,A146 ,A147 ,A14 8 , A149 ,A150 ,A151 , A15 2 25 .0 

A15 3 , A15 4 ,A15 5 ,A156 ,A15 7 ,A15 8 ,A15 9 ,A160 25 .0 

A161,A162,A163,A164,A165 ,A166 ,A167,A168 25.0 

A169 ,A170 ,A171 ,A172 ,A173 ,A174 ,A175 ,A176 25.0 

A177 ,A178 ,A179 ,A130 25.0 

A181 ,A182 , A18 3 , A184 25.0 

A185 ,A186 ,A187 ,A188 25.0 

A189 ,A190 ,A191 ,A192 25.0 



5 7 



FINAL 

25.38 
25 . 35 
25.28 
25 .21 
25 .14 
25.07 
25.00 

24.95 
24 . 89 

24.85 
24 .98 

24 .87 

24.88 

24 . 89 

24 . 89 

24.90 

24.91 
24.91 
24.91 
24.91 

24.96 
24.96 
24.96 

24.95 

24.68 

24.69 

24.71 

24.73 

24.74 

24.75 
24.74 

24.72 
24.69 
24.65 
24 . 82 

24.96 
24.96 
24 .87 

24.86 



Table XIII (cont’d) 



AREA ( sq . in . ) 



MEMBER 


INITIAL 


FINAL 


A19 3 , A19 4 , A19 5 ,A196 


25.0 


24.95 


A19 7 ,A19 8 , A19 9 ,A200 


25.0 


24.97 


A201 ,A20 2 


25.0 


24.91 


A20 3 , A204 


25 .0 


24 .92 


A205 ,A206 ,A207 ,A208 


25 .0 


24 . 86 


A209 ,A210,A211,A212 


25.0 


24 .87 


A21 3 ,A2 14 


25.0 


24 . 89 


A215 , A2 16 


25.0 


24.90 


A217 , A2 18 


25 .0 


24.91 


A219 ,A22 0 


25 .0 


24.92 


A221 ,A22 2 


25 . 0 


24.92 


A22 3 , A224 


25.0 


24 .93 


A225 ,A226 


25.0 


24.93 


A227 ,A228 


25.0 


24 .94 


A229 ,A2 30 


25 .0 


24.94 


A2 31 , A2 32 


25 .0 


24.94 


A2 3 3 , A2 34 


25.0 


24.93 







COORDINATES 


( in . ) 






JOINT 


X 


INITIAL 

Y 


Z 


X 


FINAL 

Y 


Z 


1 


120 


0* 


120 


120 .5 


0 


121.4 


5 


111 


120 


111 


110.7 


120 .4 


112.1 


9 


102 


240 


10 2 


102.2 


241.2 


10 3.2 


13 


93 


360 


93 


93.1 


362 .6 


93.7 


17 


84 


4 80 


84 


84 .0 


483.8 


84.4 


21 


75 


600 


75 


74 .9 


604.5 


75 .2 


25 


66 


720 


66 


65 .9 


724.2 


65 .9 


29 


57 


840 


57 


56 . 8 


842 .6 


56 . 8 


33 


48 


960 


48 


47 .8 


960.4 


47.8 



Initial weight = 84,524 lbs. Final weight = 83,142 lbs. 

*Final weight (AREA VARBIABLES ONLY) = 28,900 lbs. 

* This clearly indicates optimum not reached with previous case. 
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Table XIII (cont'd) 



COORDINATES (in.) 



JOINT 


X 


INITIAL 

Y 


Z 


X 


FINAL 

Y 


Z 


37 


39 


1080 


39 


38 . 8 


1077.4 


38.8 


41 


30 


0* 


30 


29 .9 


- 


29 .9 


45 


30 


1248 


30 


29 . 8 


1212 .0 


29 .9 


49 


0* 


1248 


90 


- 


1203.5 


90 .0 


51 


15 


1224 


60 


15 .0 


1194.1 


58.9 


53 


15 


1248 


60 


15 .0 


1214.7 


59 .2 



*NOT A DESIGN VARIABLE 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The finite element code was presented coupled to an 
optimizer for truss analysis and design. Trusses were 
designed for minimum weight with multiple load conditions 
considered . 

The displacement method for static analysis and the sub- 
space iteration method for eigenvalues were applied. 

Several examples were considered. In every case the code 
worked as an analysis tool, and significant weight reductions 
were obtained with the coupled optimizer CONMIN. Run times 
and results compared to the same test cases with other codes 
indiate that the code is competitive as a design tool. 

B. RECOMMENDATIONS 

The following recommendations may be of value for follow 
on work: 

1. The code should be modified so the user can access 
stresses directly. 

2. The code should be extended to include other elements 
such as frames and plates . 

3. An out of core equation solver should be added. 

4 . The method of gradient calculation should be dependent 
on specific gradients required reference 3 and reference 
4 . 
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5 . Gradients of frequency constraints would benefit from 
a more efficient algorithm reference 5 . 

6 . The need for a large scale public structural optimiza- 
tion code still exists . 
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APPENDIX A 



DATA FILES 



A. INTRODUCTION 

This appendix contains the data files used to create the 
test cases in Chapter VI. Additionally the data file for the 
user's guide in complete form is presented. 
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Table XIV 



DATA FILE 3-BAR TRUSS 



(STRESS) 



3- BAR TRUSS (TEST CASE A-STRESS ONLY) 



3 

0 

0 

1 

2 

3 

4 



4 4 

0 1 

386.4 

- 10 . 



2 

2 



0. 

0 . 

10 . 



10000000 . 
1 1 



2 

3 

1 

2 

3 

4 
1 
2 
1 



3 

4 
1 
0 
1 
1 



*1 

2 

2 

2 

1 

0 

1 

1 



1 
0 

0.0 

0. 

- 10 . 

0. 

0. 

-15000. 
1 1 



2 

1 



0 . 

0. 

0. 

0. 

20000 * 

1. 

1. 

1. 



1 

0 

0 

1 



4. 



0 

0 

0 

0 



14140. -14140. 



1. 

.01 

10 . 

10 . 

1.0 

20 . 



2 -14140. -14140. 

500. 

1. 

.01 
10 . 



0. 

0. 



0 

0 

0 

0 
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HOOH 



Table XV 



DATA FILE 3-BAR TRUSS 



( STRESS , DISPLACEMENT, BUCKLING, FREQUENCY) 



3- BAR TRUSS ( TEST CASE D-STRESS?DI SPL. t3U3 KL. tFREQ. ) 



3 

2 

0 

1 

2 

3 

4 



4 4 

2 1 

386*4 

- 10 . 



0. 

0. 

10 . 



10000000 . 



1 
2 
3 
1 
2 

3 

4 
1 
2 
1 

2 -14140 



1. 

.01 

10 . 

10 . 

1.0 

20 . 



1 

3 

4 
1 
0 
1 
1 



,1 

2 

2 

2 

1 

0 

1 

1 



1 

1 2 

0.0 

0. 

- 10 . 

0. 

0. 

-15000. 



1 

2 

1 



0. 

0. 

0. 

0. 

20000 . 

1. 

1. 

1. 



1 

0 

0 

1 



4. 



0 

0 

0 

0 



0 

0 

0 

0 



14140. -14140. 



-14140. 



0. 

0. 



20 . 



2 

2 



500. 

1 • 
.01 
10 . 



1 

2 



1 

1 



-.02 

-.02 



.02 

.02 
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HOOH 



Table XVI 



DATA FILE 18 -BAR TRUSS 
(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 



18-RAR TRUSS CAS £ (TEST 



18 
2 
0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
11 

10000000 * 



11 11 
4 1 

386.4 



1250. 

1000 . 

1000 . 

750. 

750. 

5 00. 
500. 
250. 

2 50. 

0 . 



250. 

250. 

250. 

0. 

250. 

0. 

250. 

0. 

250. 



1 

1 

0. 



CASE STRESS/ D I SP./BUC:<L IMG/ FREQ, 
1 1 0 



0 1. 



- 20000 . 



1 


1 


2 


1 


1 


10 


2 


1 


3 


1 


2 


10 


3 


2 


3 


1 


3 


10 


4 


2 


4 


1 


1 


10 


5 


3 


4 


1 


4 


10 


6 


3 


5 


1 


2 


10, 


7 


4 


5 


1 


3 


10 


8 


4 


6 


i 


1 


10 


9 


5 


6 


1 


4 


10 


10 


5 


7 


1 


2 


10 


11 


6 


7 


1 


3 


10 


12 


6 


8 


1 


1 


10 


13 


7 


8 


1 


4 


10 


14 


7 


9 


i 


2 


10 


15 


8 


9 


1 


3 


10 


16 


8 


10 


1 


1 


10 


17 


9 


10 


1 


4 


10 


18 


9 


11 


1 


2 


10 


l 


0 


0 


1 






2 






1 






3 






1 






4 






1 






5 






1 






6 






1 






7 






1 






8 






1 






9 






1 






10 


1 


1 


1 






11 


1 


1 


1 






5 












1 




— 


20000. 






2 






20000. 







4 

6 

8 

1 



5000. 



- 20000 . 

- 20000 . 

- 20000 . 



0 1 . 

1. 

1. 



1. 

1. 



1. 



20000 * 



4. 



30. 


30. 


30. 


30. 






1 . 

1000. 


1 . 

1000. ' 


1. 

1000. 


1. 

1000. 






1000. 


0. 


750. 


0. 


500 . 


0. 


700. 


-50. 


500. 


-50. 


200. 


-50. 


1200. 

1 

4.0 


200. 

2 1 


1000. 

-10. 


150. 

10. 


750. 


150. 



250, 

50. 

500, 



150 . 
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Table XVII 



DATA FILE 2 5 -BAR SPACE TOWER 
(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 



25-8AR SPACE TRUSS (TEST CASE STRESS/O I SP./ 8UCKL ING/ FREQ. ) 



25 

2 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 



10000000 * 



10 10 
2 2 
386.4 
-37.5 
37.5 
-37.5 
37.5 
37.5 
-37.5 
- 100 . 
100 . 
100 
- 100 . 



2 

2 



1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 
17 
13 

19 

20 
21 
22 

23 

24 

25 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
2 
1 
2 
4 
1 
2 
3 
6 
1 
2 

2 . 

0.00001 



1 

1 

2 

1 

2 

2 

2 

1 

1 

3 

4 

4 

5 
3 

6 
9 
8 

7 

8 

10 

9 

10 

7 

8 
9 



.1 

2 

4 

3 

5 

6 

4 

5 
3 

6 
6 

5 

3 

6 

10 

7 

4 

5 

4 
3 

5 

6 
6 

3 

4 

5 



1 
1 

0. 

0.0 

0.0 

37.5 

37.5 

-37.5 

-37.5 

100 . 

100 . 
-100. 
- 100 . 
-40000. 



1 

2 

2 

2 

2 

3 

3 

3 

3 

4 

4 

5 

5 

6 
6 
6 
6 
7 
7 
7 

7 

8 
8 
8 
8 



200 . 

200 . 

100 . 

100 . 

100 . 

100 . 

0. 

0 . 

0. 

0. 

40000. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2 . 

2 . 

2 . 

2 . 

2 . 

2 . 

2. 

2 . 

2 . 

2 . 

2 . 

2 . 

2 . 



0 

0 

1 

1 

1 

1 

4 

4 

4 

4 



0 

0 

2 

2 

2 

2 

5 

5 

- 5 

5 

39.274 



- 1 . 

1. 

1. 

- 1 . 

- 1 . 

1. 

1. 

- 1 . 



1. 

1. 

- 1 . 

- 1 . 

1. 

1. 

- 1 . 

- 1 . 



0. 

0, 

1000. 

0. 

500. 

500. 

500. 

500. 



20000 . 

- 20000 . 

10000. 
10000. 

0. 

0. 



-5000. 

-5000. 

-5000. 

-5000. 

0. 

0. 



2 . 

0.00001 



2 . 

0.00001 



2 . 

0.00001 



2 . 

0. 00001 



16.0 



50. 


50. 


50. 




50. 


50. 


37.5 


37.5 


100. 




100. 


100 


.1 


.1 


.1 




. 1 


.1 


500. 


500. 


500. 




500. 


500 


1 1 


1 -.35 




.35 






2 1 


1 -.35 




.35 






1 2 


1 -.35 




. 35 






2 2 


1 -.35 




.35 







2 . 

0.00001 

50. 



2 . 

O.OOOOX 

50. 



2 . 

0.00001 

50. 
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DATA FILE 2 34 -BAR SPACE TOWER 



234-BAR SPACE TOWER OESIGN (TEST CASE STRES S , BUCK. , OIS PL FREQ. ) 



234 


58 4 


1 5 


1 20 


0 












1 


3 2 


1 2 
















0 


386.4 


0. 
















1 


120. 


0. 


120. 


1 


0 


2 


1. 


0. 


1 . 


2 


-120. 


0. 


120. 


1 


0 


2 


-1. 


0. 


1. 


3 


-120. 


0.0 


-120. 


1 


◦ 


2 


-1. 


0. 


-1. 


4 


120. 


0. 


-120. 


1 


0 


2 


u 


0. 


-1. 


5 


111. 


120. 


111. 


3 


4 


5 


i. 


1. 


1. 


6 


-111. 


120. 


111. 


3 


4 


5 


-l. 


1. 


1 . 


7 


-111. 


120. 


-111. 


3 


4 


5 


-i. 


1. 


-1. 


8 


111. 


120. 


-111. 


3 


4 


5 


i. 


1. 


-1. 


9 


102. 


240. 


102. 


6 


7 


8 


i. 


1. 


1 . 


10 


-102. 


240. 


102. 


6 


7 


8 


-i. 


1. 


1. 


11 


-102. 


240. 


-102. 


6 


7 


8 


-i. 


1. 


-1. 


12 


102. 


240. 


-102. 


6 


7 


8 


i. 


1. 


-1. 


13 


93 . 


360. 


93. 


9 


10 


11 


l. 


1. 


1. 


14 


-93. 


360. 


93. 


9 


10 


11 


-i. 


1. 


1. 


15 


-93. 


360. 


-93. 


9 


10 


11 


-l. 


1. 


-1. 


16 


93. 


360. 


-93. 


9 


10 


11 


i. 


1. 


-1. 


17 


84. 


480. 


84. 


12 


13 


14 


i. 


1. 


1. 


18 


-84. 


480. 


84. 


12 


13 


14 


-l. 


1. 


1. 


19 


-84. 


480. 


-84. 


12 


13 


14 


-l. 


1. 


-1. 


20 


84. 


480. 


-84. 


12 


13 


14 


i. 


1. 


-1. 


21 


75. 


600. 


75. 


15 


16 


17 


i. 


1. 


1. 


22 


-75. 


600. 


75. 


15 


16 


17 


-i. 


1. 


1. ' 


23 


-75. 


600. 


-75. 


15 


lo 


17 


-l. 


1. 


— 1 • 


24 


75. 


600. 


-75. 


15 


16 


17 


i. 


1. 


-1 . 


25 


66. 


720. 


66. 


18 


19 


20 


i . 


1. 


1. 


26 


-66. 


720. 


66. 


18 


19 


20 


-i. 


1. 


1 . 


27 


-66. 


720. 


-66. 


18 


19 


20 


-i. 


1. 


-1. 


28 


66. 


720. 


-66. 


16 


19 


20 


i. 


1. 


. -1. 


29 


57. 


840. 


57. 


21 


22 


23 


i. 


1. 


1. 


30 


-57. 


840. 


57. 


21 


22 


23 


-l. 


1. 


1. 


31 


-57. 


840. 


-57. 


21 


22 


23 


-l. 


1. 


-1. 


32 


57. 


640. 


-57. 


21 


22 


23 


i. 


1. 


-1. 


33 


48. 


960. 


48. 


24 


25 


26 


i. 


1. 


1. 


34 


-48. 


960. 


48. 


24 


25 


26 


-i. 


1. 


1. 


35 


-48. 


960. 


-46. 


24 


25 


26 


-l. 


1. 


-1. 


36 


48. 


960. 


-48. 


24 


25 


26 


i. 


1. 


-1. 


37 


39. 


1080. 


39. 


27 


28 


29 


l. 


1. 


1. 


38 


-39. 


1080. 


39. . 


27 


28 


29 


-i. 


1- 


1. 


39 


-39. 


1080. 


-39. 


27 


28 


29 


-i. 


1. 


-1. 


40 


39. 


1080. 


-39. 


27 


28 


29 


l. 


1. 


-1. 


41 


30. 


1200. 


30. 


39 


0 


40 


i. 


0 


1. 


42 


-30. 


1200. 


30. 


39 


0 


40 


-i. 


0 


1. 


43 


-30. 


1200. 


-30. 


39 


0 


40 


-l. 


0 


-1. 


44 


30. 


1200. 


-30. 


39 


0 


40 


i. 


0 


-1. 


45 


30. 


1248. 


30. 


41 


30 


42 


i. 


1. 


1. 


46 


-30. 


1248. 


30. 


41 


30 


42 


-i. 


1. 


1. 


47 


-30. 


1248. 


-30. 


41 


30 


42 


-i. 


1. 


-1. 


48 


30. 


1248. 


-30. 


41 


30 


42 


i. 


1. 


-1. 


49 


90. 


1248. 


0. 


0 


31 


32 


i. 


1. 


0 


50 


-90. 


1248. 


0. 


0 


31 


32 


-i. 


1. 


0 


51 


60. 


1224. 


15. 


33 


34 


35 


i. 


1. 


1. 


52 


-60. 


1224. 


15. 


33 


34 


35 


-i. 


1* 


1. 


53 


60. 


1248. 


15. 


36 


37 


38 


i. 


1. 


1. 


54 


-60. 


1248. 


15. 


36 


37 


38 


-i. 


1. 


1. 


55 


-60. 


1224. 


-15. 


33 


34 


35 


-l. 


1. 


-1. 


56 


60. 


1224. 


, -15. 


33 


34 


35 


l. 


1. 


-1. 


57 


-60. 


1248. 


; -15. 


36 


37 


38 


-i. 


1. 


-1. 


58 


60. 


1248. 


• -15. 


36 


37 


38 


i. 


1. 


-1. 


30000000. .1 


-15000, 


. 20000. 




4. 










1 


1 5 


1 I 


10. 














2 


2 6 


1 1 


10. 














3 


3 7 


1 1 


10. 














4 


4 8 


1 1 


10. 














5 


5 9 


1 2 . 


10. 














6 


6 10 


1 2 


10. 














7 


7 11 


1 2 


10. 














8 


8 12 


1 2 


10. 














9 


9 13 


1 3 


10. 














10 


10 14 


1 3 


10. 






V 






- 


11 


11 15 


1 3 


10. 
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Table XVIII (cont’d) 



12 


12 


16 


1 


3 


10 . 


13 


13 


17 


1 


4 


10 . 


14 


14 


13 


1 


4 


10 . 


15 


15 


19 


1 


4 


10 . 


16 


16 


20 


1 


4 


10 . 


17 


17 


21 


1 


5 


10 . 


18 


18 


22 


1 


5 


1 0 . 


19 


19 


23 


1 


5 


10 . 


20 


20 


24 


1 


5 


10 . 


21 


21 


25 


1 


6 


10 . 


22 


22 


26 


1 


6 


10 . 


23 


23 


27 


1 


6 


10 . 


24 


24 


28 


1 


6 


10 . 


25 


25 


29 


1 


7 


10 . 


26 


26 


30 


1 


7 


10 . 


27 


27 


31 


1 


7 


10 . 


28 


28 


32 


1 


7 


10 . 


29 


29 


33 


1 


8 


10 . 


30 


30 


34 


1 


8 


10 . 


31 


31 


35 


1 


8 


10 . 


32 


32 


36 


1 


8 


10 . 


33 


33 


37 


1 


9 


10 . 


34 


34 


38 


1 


9 


10 . 


35 


35 


39 


1 


9 


10 . 


36 


36 


40 


1 


9 


10 . 


37 


37 


41 


1 


10 


10 . 


38 


38 


42 


1 


10 


10 . 


39 


39 


43 


1 


10 


10 . 


40 


40 


44 


1 


10 


10 . 


41 


41 


45 


1 


11 


10 . 


42 


42 


46 


1 


11 


10 . 


43 


43 


47 


1 


11 


10 . 


44 


44 


48 


1 


11 


10 . 


45 


5 


6 


1 


12 


10 . 


46 


6 


7 


1 


12 


10 . 


47 


7 


8 


1 


12 


10 . 


48 


8 


5 


1 


12 


10 . 


49 


9 


10 


1 


13 


10 . 


50 


10 


11 


1 


13 


10 . 


51 


11 


12 


1 


13 


10 . 


52 


9 


12 


1 


13 


10 . 


53 


13 


14 


1 


14 


10 . 


54 


14 


15 


1 


14 


10 . 


55 


15 


16 


1 


14 


10 . 


56 


13 


16 


1 


14 


10 . 


57 


17 


18 


1 


15 


10 . 


58 


18 


19 


1 


15 


10 . 


59 


19 


20 


1 


15 


10 . 


60 


17 


20 


1 


15 


10 . 


61 


21 


22 


1 


16 


10 . 


62 


22 


23 


1 


16 


10 . 


63 


23 


24 


1 


16 


10 . 


64 


21 


24 


1 


16 


10 . 


65 


25 


26 


1 


17 


10 . 


66 


26 


27 


1 


17 


10 . 


67 


27 


28 


1 


17 


10 . 


68 


25 


28 


1 


17 


10 . 


69 


29 


30 


1 


18 


10 . 


70 


30 


31 


1 


18 


10 . 


71 


31 


32 


1 


18 


10 . 


72 


29 


• 32 


1 


18 


/ 10 . 
/ 10 . 


73 


33 


34 


1 


19 


74 


34 


35 


1 


19 


10 . 


75 


35 


36 


1 


19 


10 . 


76 


33 


36 


1 


19 


10 . 


77 


37 


38 


1 


20 


10 . 


78 


33 


39 


1 


20 


10 . 


79 


39 


40 


1 


20 


10 . 


80 


37 


40 


1 


20 


10 . 


81 


41 


42 


1 


21 


10 . 


82 


42 


43 


1 


22 


10 . 


83 


43 


44 


1 


21 


10 . 


84 


41 


44 


1 


22 


10 . 


85 


45 


46 


1 


23 


10 . 
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86 


46 


47 


1 


24 


10 


87 


47 


48 


1 


23 


10, 


88 


45 


48 


1 


24 


10 


89 


1 


6 


1 


2 5 


10 


90 


2 


5 


1 


25 


10 


91 


2 


7 


1 


25 


10 


92 


3 


6 


1 


25 


1 0 , 


93 


3 


8 


1 


2 5 


10, 


94 


4 


7 


1 


25 


10 


95 


4 


5 


1 


25 


10, 


96 


1 


8 


1 


25 


10 


97 


5 


10 


1 


26 


1 0 


98 


6 


9 


1 


26 


10, 


99 


6 


11 


1 


26 


10 


100 


7 


10 


1 


26 


10 


101 


7 


12 


1 


26 


10 


102 


8 


11 


1 


26 


10 


103 


8 


9 


1 


26 


10 


104 


5 


12 


1 


26 


10, 


105 


9 


14 


1 


27 


10 


106 


10 


13 


1 


27 


10, 


107 


10 


15 


1 


27 


10, 


108 


11 


14 


1 


27 


10 


109 


11 


16 


1 


27 


1 0 , 


no 


12 


15 


1 


27 


10, 


in 


12 


13 


1 


27 


10 


112 


9 


16 


1 


27 


10, 


113 


13 


18 


1 


28 


10 


114 


14 


17 


1 


28 


10 


115 


14 


19 


1 


28 


10, 


116 


15 


18 


1 


28 


10, 


117 


15 


20 


1 


28 


10, 


118 


16 


19 


1 


28 


10, 


119 


16 


17 


1 


28 


10 


120 


13 


20 


1 


28 


10. 


121 


17 


22 


1 


29 


10, 


122 


18 


21 


1 


29 


10 


123 


18 


23 


1 


29 


10, 


124 


19 


22 


1 


29 


10 


125 


19 


24 


1 


29 


10, 


126 


20 


23 


1 


29 


10, 


127 


20 


21 


1 


29 


10, 


128 


17 


24 


1 


29 


10, 


129 


21 


26 


1 


30 


10 


130 


22 


25 


1 


30 


10 


131 


22 


27 


1 


30 


10, 


132 


23 


26 


1 


30 


10 


133 


23 


28 


1 


30 


10 


134 


24 


27 


1 


30 


10 


135 


24 


25 


1 


30 


10 


136 


21 


28 


1 


30 


10, 


137 


25 


30 


1 


31 


10 


138 


26 


29 


1 


31 


10 


139 


26 


31 


1 


31 


10 


140 


27 


30 


1 


31 


10 


141 


27 


32 


1 


31 


1 0 , 


142 


28 


31 


1 


31 


10 


143 


28 


29 


1 


31 


10 


144 


25 


32 


1 


31 


10 


145 


29 


34 


1 


32 


10 


146 


30 


33 


1 


32 


1 0 


147 


30 


35 


1 


32 


1 0 


148 


31 


34 


1 


32 


10 


149 


31 


36 


1 


32 


10 


150 


32 


35 


1 


32 


10 


151 


32 


33 


1 


32 


10 


152 


29 


36 


1 


32 


10 


153 


33 


38 


1 


33 


10 


154 


34 


37 


1 


33 


10 


155 


34 


39 


1 


33 


10 


156 


35 


38 


1 


33 


10 


157 


35 


40 


1 


33 


1 0 


158 


36 


39 


1 


33 


10 


159 


36 


37 


1 


33 


10 
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33 

37 

39 

38 

39 

39 

40 

40 

37 

41 

42 

42 

43 

43 

44 

44 

41 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

50 

49 

49 

50 

50 

49 

45 

46 

47 

48 

51 

52 

55 

56 

53 

54 

51 

52 

53 

47 

45 

46 

43 

41 

42 

44 

5 

6 

9 

10 

13 

14 

17 

13 

21 

22 

25 

26 

29 

30 

33 

34 

37 

38 

41 

42 

45 



( cont r d) 



40 


1 


33 


10. 


42 


1 


34 


10. 


41 


1 


34 


10. 


43 


1 


34 


10. 


42 


1 


34 


10. 


44 


1 


34 


10. 


43 


1 


34 


10. 


41 


1 


34 


1 0 . 


44 


1 


34 


10. 


46 


1 


35 


10. 


45 


1 


35 


10. 


47 


1 


35 


10. 


46 


1 


35 


10. 


48 


1 


35 


10. 


47 


1 


35 


10. 


45 


1 


35 


10. 


48 


1 


35 


10. 


51 


1 


36 


10 . 


52 


1 


36 


10. 


55 


1 


36 


10. 


56 


1 


36 


10. 


53 


1 


37 


10. 


54 


1 


37 


10. 


57 


1 


37 


10. 


58 


1 


37 


10. 


51 


1 


38 


10. 


52 


1 


38 


10. 


55 


1 


38 


10. 


56 


1 


38 


10. 


53 


1 


39 


10. 


54 


1 


39 


10. 


57 


1 


39 


10. 


58 


1 


39 


10. 


51 


1 


40 


10. 


52 


1 


40 


10. 


55 


1 


40 


10. 


56 


1 


40 


10. 


53 


1 


41 


10. 


54 


1 


41 


10. 


57 


1 


41 


10. 


58 


1 


41 


10. 


58 


1 


42 


10. 


57 


1 


42 


10. 


56 


1 


43 


10. 


55 


1 


43 


10. 


48 


1 


44 


10. 


54 


1 


44 


10. 


58 


1 


44 


10. 


57 


1 


44 


10. 


52 


1 


45 


10. 


5o 


1 


45 


10. 


55 


1 


45 


10. 


51 


1 


45 


10. 


7 


1 


46 


10. 


8 


1 


46 


10. 


11 


1 


47 


10. 


12 


1 


47 


10. 


15 


1 


48 


10. 


16 


1 


48 


10. 


19 


1 


49 


10. 


20 


1 


49 


10. 


23 


1 


50 


10. 


24 


1 


50 


10. 


27 


1 


51 


10. 


28 


1 


51 


10. 


31 


1 


52 


10. 


32 


1 


52 


10. 


34 


1 


53 


10. 


36 


1 


53 


1 0 . 


39 


1 


54 


10. 


40 


1 


54 


10. 


43 


1 


55 


10. 


44 


1 


55 


10. 


47 


1 


56 


10. 
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234 

1 

2 

3 

4 
2 



46 

1 

1 

1 

1 



48 

1 

1 

1 

1 



56 



10 , 



49 


6000. 


-20000. 


0.0 










50 


6000. 


-20000. 


0.0 










Z 

49 


6000. 


-20000. 


0.0 










50 


-6000. 


-20000. 


0.0 










Z 

49 


6000. 


-20000. 


0.0 


* 








50 

9 


3000. 


-10000. 


5 000. 










z 

49 


3000. 


-10000. 


-5000. 










50 

9 


3000. 


-10000. 


5000. 










z 

49 


-3000. 


10000. 


5000. 










50 


-3000. 


10000. 


-5000. 






1 




49 


200. 














50 


200. 














;. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


i. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


i. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


>. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


>. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


>. 


25. 


25. 


25. 


25. 


25. 


25. 


25 


>. 


25. 


25. 


25. 


25. 


25. 


25. 


25 



0*000001 .000001 
0.000001 .000001 



0.000001 

0.000001 

.000001 

.000001 

.000001 



.000001 

.000001 

.000001 

.000001 

.000001 



.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 



.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 



.000001 

.000001 

.000001 

.000001 

.000001 

.000001 

.000001 



.000001 
.000001 
.000001 
.000001 
.000001 
►000001 
►000001 



. 000001 
. 000001 
.000001 
. 000001 
.000001 
.000001 
.000001 



100. 


100. 


100. 


100. 


100. 


100. 


100. 


100 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100, 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


ICO, 


100. 


100. 


100. 


100.- 


100. 


100. 


100. 


100 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100 


120. 


120. 


111. 


120. 


111. 


102. 


240. 


102 


93. 


360. 


93. 


84. 


480. 


84. 


75 . 


600 


75. 


66. 


720. 


66 • 


57. 


840. 


57. 


48. 


960. 


48. 


39. 


1080. 


39 • 


1248. 


1248. 


90. 


15. 


1224. 


60. 


15. 


1248. 


60. 


30. 


30. 


30. 

0. 


30. 

0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 



0. 



0. 



.000001 
. 0000C1 
.000001 
.000001 
.000001 
.000001 
.000001 



2000. 


2000. 


2000. 


2000. 


2030. 


2000. 


2000. 


2000 


2000. 


2000. 


2000. 


2000. 


2000. 


2000. 


2000. 


2000 


2000. 


2000. 


2000. 


2000. 


2000. 


2C00 • 


* 2000. 


2000 


2000. 


2000. 


2000. 


2000. 


2000. 


2000. 


2000. 


2000 


2000. 


2000. 


2000. 


, 2000. 


2000. 


2000. 


2000. 


2000 


2000. 


2000. 




/ 










49 


1 


1 


-.3 


.3 










50 


1 


1 


-.3 


.3 










49 


1 


2 


-.3 


.3 










50 . 


1 


2 


-.3 


.3 










49 


1 


3 


-.3 


.3 










50 


l 


3 


-.3 


. 3 










49 


1 


4 


-.3 


.3 










50 


1 


4 


-.3 


.3 










49 


1 


5 


-.3 


.3 










50 


1 


5 


-.3 


.3 










49 


2 


1 


-.3 


.3 




\ ' 






50 


2 


1 


-.3 


.3 






- 


49 


2 


2 


-.3 


.3 
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Table XVIII (ccnt T d) 



4,5 



50 

49 

50 

49 

50 

49 

50 



2 

2 

2 

2 

2 

2 

2 



2 

3 

3 

4 

4 

5 
5 



-.3 

-.3 

-•3 

-.3 

-.3 

-.3 

-.3 
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APPENDIX B 



PROGRAM ORGANIZATION 



A. DESCRIPTION 

The program organization is laid out in the following 
flow charts . The main driver program is arranged to call a 
subdriver, XMSADT , and the optimizer of the user's choice. 
All changes for replacing an optimizer occur in MSADT . This 
allows for easy testing of several optimizers on the same 
problem. 

XMSADT may be called from the main for input, analysis, 
and output . Printed output may vary as the user requires . 

A complete listing of all subroutines and their functions is 
given in Table 
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Table XIX 



ANALYSIS CONTROL 



ANALYSIS ICALC = 2 



ICALC=2 
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Table XX 



OUTPUT CONTROL 



OUTPUT FOR INFOsO/1/2 



ICALC= 3 




(PRINT) 
-GRADIENTS 
SPECIFIED BY 
INPUT CONTROL 
(IDVCLC) 
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Table XXI 

PROGRAM FLOW CHART 



MSADT 




OPTIMIZER 



(INPUT) 
SAD TO 1 



SADT46 

SADT02 



SADT11 



SADTO 4- 

SADT07 

SADT08 



SADT03 

SADT05 



SADT61 - 



SADT50 

SADT49 



— SADT06 
-SADT43 
— SADTO 5 



SADT1| 



SADTO 8 
• EIGRQ1 



SADT1-5 — SADT16 
■SADT2 3 

canalysTs) 



EIGR02 EIGR03 



SADT7 3- 

SADT75 

SADT76 

SADT77 



SADT36 
SADT2 0“ 

SADT21 



SADTO 3 
SADT4 8 

SADT37 



(GRADIENTS) 



4 =- 



SADT74 

SADT53 

SADTQ9 

SADT7 9- 
SADT17 
SADT62 

SADT71 



(OUTPUT) 
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TABLE XXII 



FILENAME 

EIGR01 
EIGR02 
EIGR03 
SADTO 1 

SADT02 

SADT03 

SADTO 4 

SADTO 5 

SADT06 

SADT07 

SADT08 

SADT09 

SADTIO 
SADT1 1 
SADT14 

SA DTI 5 

SADT16 
SA DT 17 
S A DTI 8 
SADT19 



SUBROUTINE USES 



TITLE 



SOLVES EIGENVALUE PROBLEM | A-ALAMBD ft B | X=0 

SOLVES EIGENVALUE PROBLEM 
SOLVES EIGENVALUE PROBLEM 

THIS ROUTINE READS AND PRINTS INPUT DATA 
AND ORGANIZES PSEUDO-DYNAMIC STORAGE 
ALLOCA TION 

BUILDS VECTORS JC AND UK FOR FINITE ELEMENT 
STRUCTURAL ANALYSIS 

BUILDS THE ELEMENT STIFFNESS MATRIX FOR 
TRUSS FINITE ELEMENT NO. II 

BUILDS ELEMENT STIFFNESS MATRICES AND ADDS 
THEM TO THE GLOBAL STIFFNESS MATRIX 

SUPERIMPOSES THE ELEMENT STIFFNESS VECTOR EK 
OF ELEMENT 1 1 ON THE GLOBAL VECTOR AK 

BUILDS GLOBAL LUMPED MASS MATRIX 

LU DECOMPOSES SYMMETRIC, POSITIVE-DEFINITE, 
SPARCE MATRICES, THE UPPER TRIANGLE OF WHICH 
IS STORED IN VECTOR AK WITH LEADING ZEROES 
OMITTED 

FORWARD AND BACK SUBSTITUTES TO YIELD A SOLN 
TO A SET OF LINEAR EQNS ( DECOMPOSED W/ 

SADT07) 

PRINTS ALL JOINT DISPLACEMENTS FOR EACH 
LOAD CONDITION OF A FINITE ELEM. STRUCTURE 

PERFORMS TRUSS FIXED-GEOMETRY DESIGN 

ROUTINE TO ORGANIZE ANALYSIS 

CALCULATES PARTIAL DERIVATIVE OF TRUSS ELEM 
K MATRIX WRT COORD DOF II AT NODE J1 

CALCULATESS VALUES OF ALL DESIGN AND BEHAV- 
IORIAL CONSTRAINTS AS DEFINED BY PROGRAM 
"SAD" 

CALCULATES STRESS IN BAR ELEMENT II UNDER 
LOAD CONDITION JJ 

PRINTS FORCES OR STRESSES IN BAR ELEMENTS 



CALCULATES GRADIENT INFO IN FINITE ELEM 
STRUCTURAL ANALYSIS 

ADDS ELEMEMENT MASS AA OF ELEMEMENT II 
TO GLOBAL MASS MATRIX 
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SADT20 


CALCULATES THE PARTIAL DERIVATIVES OF THE 
DISPLACEMENTS WRT INDEPENDENT DESIGN 
VARIABLE II 


SADT21 


CALCULATES STRESS AND GRADIENT OF STRESS 
IN BAR ELEM II UNDER LOAD COND JJ GIVEN 
DISPLACEMENT VECTOR PL AND GRADIENTS OF 
THE DISPLACEMENTS IN VECTOR DU 


SADT22 


PERFORMS STRESS-RATIO DESIGN OF A FIXED 
GEOMETRY TRUSS 


SADT23 


CALCULATES HEIGHT OF TRUSS GIVEN AREAS AND 
LENGTHS 


SADT26 


CALCULATES THE GRADIENT OF WEIGHT WRT 
RECIPROCAL DESIGN VARIABLES 


SADT31 


CALCULATES GRADIENT INFORMATION AND 
STORES IN ARRAY TAY 


SADT34 


CONMIN EXTERNAL FOR TRUSS FIXED GEOMETRY 
DESIGN 


SADT35 


CALCULATES GRADIENT INFORMATION AND 
STORES IN ARRAY TAY 


SADT36 


CALCULATES ( XEIG-T* AM* XEIG) FOR 
GRADIENT CALCULATIONS IN FREQUENCY 
CONSTRAINTS 


SADT37 


CALCULATES EIGENVALUE GRADIENT INFORMATION 
IN FINITE ELEMENT STRUCTURAL ANALYSIS AND 
DESIGN 


SADT40 


CALCULATES GRADIENTS OF DISPLACEMENTS 
WRT COORDINATE DESIGN VARIABLE II 


SADT41 


CALCULATES GRADIENT INFORMATION WRT 
COORDINATE VARIABLES 


SADT42 


CALCULATES GRADIENT OF WEIGHT WRT 
COORDINATE DESIGN VARIABLES 


SADT43 


CALCULATES GRADIENTS OF OBJECTIVE AND 
ACTIVE CONSTRAINTS IN OPTIMUM GEOMETRY 
DIRECTION FINDING PROBLEM 


SADT44 


PERFORMS TRUSS FIXED GEOMETRY DESIGN AND 
EVALUATES CONSTRAINTS 


SADT45 


PERFORMS TRUSS OPTIMUM GEOMETRY DESIGN 


SADT46 


READS INPUT INFORMATION FOR BAR ELEMENTS 


SADT48 


BUILDS THE ELEMENT MASS MATRIX FOR A 
TRUSS FINITE ELEMENT 


SADT49 


SOLVES REAL EIGENVALUE PROBLEMS USING 
THE SUBSPACE ITERATION METHOD 


SADT50 


BUILDS INITIAL SET OF BASIS VECTORS 
FOR El GEN SOL UTION BY REDUCED BASIS 
METHOD 


SADT51 


CALCULATES STRESS AND GRADIENT OF STRESS 
IN BAR ELEMENT II UNDER LOAD CONDITION 
JJ, GIVEN DI SLPL ACEMENTS CONTAINED IN 
VECTOR PL AND GRADIENT OF DISPLACEMENTS 
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IN VECTOR DU 



SADT52 

SADT53 

SADT54 

SADT6 1 
SADT62 

SADT7 1 
SADT73 
SADT74 
SADT75 

SADT76 

SADT77 
SADT79 
MS ADT 

XMSADT 



PERFORMS 1-DIMENSIONAL SEARCH IN 
COORDINATE DESIGN SPACE 

PRINTS MEMBER INFORMATION FOR BAR ELEM 

CALCULATES AND PRINTS CENTER-OF-GRA VITY 
AND INERTIAL PROPERTIES OF FINITE ELEM 
STRUCTURE 

BUILDS GLOBAL MASS MATRIX 

PRINTS NEIG EIGENVALUES STORED IN 
EIGVAL, AND THEIR CORRESPONDING 
EIGENVECTORS STORED IN XEIG 

PRINTS G VECTOR OF CONSTRAINTS 

CALCULATES GRADIENT INFORMATION 

PRINTS TAY ARRAY OF GRADIENTS 

CALCULATES GRADIENT OF WEIGHT WITH 
RESPECT TO AREA DESIGN VARIABLES 

CALCULATES GRADIENT OF WEIGHT WITH 
RESPECT TO COORDINATE DESIGN VARIABLES 

CALCULATES GRADIENT INFORMATION 

PRINTS COORDINATE INFORMATION 

DRIVER PROGRAM FOR USING THE ABOVE 
SUBROUTINES. MSADT MAY BE COUPLED 
TO OPTIMIZER OF USER'S CHOICE. 

SUBDRIVER PROGRAM FOR COUPLING 
SADT ROUTINES TO MSADT 
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